Raman Scattering cross section of Spin Ladders 
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The Raman scattering spectra from magnetic excitations in an antiferromagnetic spin-i two leg 
ladder is investigated for weak and strong interladder coupling. In the first case, a cusp in the 
Raman intensity is obtained at a frequency twice the gap. In the second case, a peak at twice the 
gap replaces the cusp. We discuss the relevance of our calculation to recent experiments on CaV^Os 
and Sri4Cu2404i- 
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Raman scattering is an experimental technique that 
has provided valuable informations about the spin dy- 
namics in quasi-one dimensional antiferromagnets in the 
recent years. Since Raman scattering is sensitive to sin- 
glet excitations, this technique is complementary of neu- 
tron diffraction which is sensitive to triplet excitations. 
It has been used to ffioie spin 1/2 chainsEnSRin 1 chainsa, 
spin Peierls systemsETEl, and spin laddersErEI. In particu- 
lar, Raman scattering has been very useful in the analy- 
sis of magnetic excitations in the spin-Peierls compound 
CuGeC>3. The bosonized theory of dimerized spin 1/2 
states that is believed to describe the dimerized low tem- 
perature phase of spin-Peierls systems predicts the ap- 
pearance of a singlet bound state of two triplet-excita- 
tions at an energy \/3A, where A is the spin gaplii Such 
singlet bound state has been successfully observed in Ra- 
man scattering experiments on CuGeC>3 for T < Tgp 
at an energy 1.79A, close to the theoretical prediction. 
Moreover, the peak was not observed in the uniform 
phase(T > Tgp) n| showing that it is characteristic of the 
dimerized phaseoO. For T > n Tgp, a broad band of mag- 
netic excitations is observedo. The theoretical analysis 
of magnetic Raman |-scattering is based on the Fleury- 
Loudon Hamiltonianl3t3, that describes the interaction 
of photons with magnetic excitations. There exists at 
present a certain amount of litterature on the theory 
of Raman— s catt ering from dimerized spin chains, both 
analyticallD'ESlla and numerically. The_case of frustrated 
spin chains has also been investigatedEZI, in relation with 
the Raman spectra of CuGeC>3 at T > Tgp. The theory 
of Ref. [I?] reproduced well the features of the spectrum 
at T > TgpQ. The application of Raman scattering to 
probej-the singlet excitations of two leg ladders is more 
recenlOij. In spin ladder systems, magnetic peaks in the 
Raman intensity were observed at twice the spin gap. 
From the theoretical point of view, some numerical cal- 
culations are availablets, but no analytic expression of 
the Raman intensity has been derived so far. In order to 
fill this gap, we discuss in the present work the Raman 



spectrum of an antiferromagnetic spin-1/2 ladder. After 
recalling some basic results on the Fleury-Loudon the- 
ory of magnetic Raman scattering, we will consider first 
the Majorana fermion approach valid for weak coupling 
and then the Bond Operator Technique (BOT) valid for 
the strong coupling case. The Majorana fermions ap- 
proach leads to a cusp in the Raman intensity at twice 
the gap, in disagreement with experiment. We discuss 
briefly what could be missing in the Majorana fermions 
description. On the other hand, the BOT predicts cor- 
rectly the presence of peaks in the Raman intensity at 
twice the gap. 

We consider two coupled antiferromagnetic S — 1/2 
Heisenberg chains, whose Hamiltonian is 
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where J > and J± > denotes the intra- and 
inter- chain antiferromagnetic interactions, respectively. 
The interaction of light with the antiferromagnetic fluc- 
tuations is described by Loudon-Fleury'al3il3 photon- 
induced super-exchange operator 
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where E i(E g) are the incident (scattered) electric field, 
and 6 y is a unit vector connecting the sites i and j, at 
which the spins "S^and S j are located. A derivation of 
(^) stapling from the Hubbard Hamiltonian can be found 
in RefH 

The Raman cross sectionEjO can be expressed as a 
function of the retardated Raman response function as: 
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ui\ and u>2 are the frequencies of the incoming and scat- 
tered radiation, respectively, u> = u>2 — Wij n\ and are 
the respective refractive index. V is the volume of the 
crystal and c the velocity of light. The retardated linear 
response function is defined as: 

/>oo 

xT M = J / e^+^'Tr {Z-^-^ ff„(0)]} , 
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where Z = Tre~^ ff and H R is the Loudon-Fleury Hamil- 
tonian (|J). 

By inserting the resolution of identity in (Q) , the Ra- 
man intensity can be written as 
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where | ^ n (m)) are eigenstates with energies E n i m \- Such 
formula can be easily interpreted as a Fermi golden rule 
averaged over the Boltzmann weight. To get informations 
on two-magnons scattering processes we should perform 
a symmetry analysis of the matrix elements appearing 
in(^), and discuss selection rules. Since the spin ladder 
Hamiltonian is invariant under translation along the legs, 
SU(2) rotation, and mirror along the leg direction, an 
eigenstate should be characterized by a (lattice) momen- 
tum defined modulo 2tt /a (where a is the lattice spacing), 
a spin and its parity under leg exchange. The Raman 
operator defined in (|2|) is rotationally and translationally 
invariant, and still invariant under leg exchange. As a 
result, the selection rules impose that the states | ^f n ) 
and | ty m ) have the same spin, momentum and parity 
under leg exchange. This implies, in particular, that at 
T = 0, transitions will only take place to states of total 
momentum zero, spin zero and same parity as the ground 
state. Let us now turn to concrete calculations. We con- 
sider the scattering for Ej and Es, parallel to the rung 
direction, thus we have 
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In the following, we will evaluate the Raman intensity 
in the weak coupling and in the strong .-coupling limit 
using the standard Matsubara techniquet3 to calculate 
the correlator xr(u). 

To evaluate the time ordered Raman susceptibility for 
the weakly coupled chains, we will employ the Majo- 
rana fcrmion representation of the spin-ladder Hamilto- 
nian (|lf) introduced by Shelton, Nersesyan and Tsvelik in 
Ref. E3l The effective Hamiltonian is expressed in terms 
of four interacting Majorana fermions. They comprise 
a degenerate triplet Q(x) (v = Left, i?ight) with bare 
mass m t = m = J± and a singlet, p v {x) with bare mass 
m s = —'dm. It has been argued in Ref. [2^ that the ef- 
fect of interactions was merely to renormalize the bare 



masses, so that interactions could be neglected. With 
this approximation, the spin ladder is described by the 
following effective Hamiltonian: 
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where p stands for the triplet or singlet mass, and C is the 
corresponding triplet or singlet operator. The thermal 
Green's function for the left and right moving triplet and 
singlet Majorana fermions are defined by: 



whose explicit expressions are 
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G RR {k, iuj n ) = G1 L {-k, iu n ) = 
G RL (k,iu>„) = Gl R {k,iuj n )* = 
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where a stands for t (triplet) or s(singlet), and uj n = 
(2n + l)n/(3 are the fermion Matsubara frequencies. 
In terms of Majorana fermions, the Raman operator 
7^7Jj Si,j • ~S 2,i, with 7 a constant, is expressed by 



Hr = Jt£R€L + IsPRpL, 
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where j t — m t j and 7 S = m s j. To arrive at this expres- 
sion, the marginal term, already neglected in the deriva- 
tion of the Hamiltonian (0), has been discarded. Injecting 
this expression into the definition of the Raman suscepti- 
bility (Q) and applying Wick's theorem, the time ordered 
expectation value at finite temperature can be written 
as: 
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-G RR (q,iu n )Gl L {-q,i{w n -v n ))}. (11) 

Explicitly, we have to compute the following integral 
and sum over the Matsubara frequencies: 
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In order to evaluate the Matsubara sum in (12), we 
have to determine the residues of the four poles of the 
expression ( |l2] ) and multiply every residue with the value 
of the Fermi function uf{z) — l/(exp(/3z) + l) at the pole. 
Adding the four terms together yields 
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where we have introduced the notation e a (q) = 
\J (vq) 2 + m, 2 . Thus performing the analytic continua- 
tion (iu> n — > lij + i0+), taking the imaginary part and 
performing the integral over q, we finally get: 



IrnxflH = 7T tanh ( 
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Formula (|1J) implies the existence of a cusp singular- 
ity in the Raman intensity at twice the spin gap due 
to the triplet excitations and another singularity at six 
times the spin gap due to the singlet modes. As re- 
sult, the noninteracting Majorana fermions representa- 
tion does-jiot reproduce the Raman peak experimentally 
observedQcl. The spectra predicted by (Q) is plotted 
on figure |l|. The absence of signal for to smaller than 
twice the gap is, in qualitative agreement with numeri- 
cal simulationsEj. It would be interesting to determine 
whether treating properly the interactions between the 
Majorana Fermions can reproduce the experimental peak 
at twice the spin gap. We now turn to a strong-coupling 
analysis of the Ramaa-susceptibility, using the bond op- 
erator representation^ ofj-auantum S=l/2 spins used by 
Gopalan, Rice and SigristEa in their mean field approach 
to spin ladders. In this representation, one starts from 
weakly coupled rungs and introduces on each rung a sin- 
glet and three triplets t' a (a — x, y, z) boson creation 
operators, that span the Hilbert space of a single rung 
when acting on a vacuum state. Since the rung can be 
in either the singlet or one of the triplet states, the con- 
dition : 
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has to be satisfied by the physical states. The repre- 
sentation of the spins Si and S2 in terms of— these sin- 
glet and triplet operators, is derived in Rcf.EilO. Sub- 
stituting this operator representation of spins into the 
original Hamiltonian, one ends up with an Hamiltonian 
quartic in boson fields. Treating the singlet operator in 
a mean field approximation and neglecting interactions 
among the triplets, one obtains the following Hamilto- 
nian quadratic in triplet operatorsci 
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The chemical potential term fi guarantees that the cond- 
tion (p5J) is satisfied on average. This Hamiltonian can be 



solved by Green's function method. One, first, introduces 
the four Green's functions Gi_ a (r) = —{T T ti ta (r)t\ Q (0)), 

-<T r f, ia (r)to, a (0)>, F} a {r) = - (T T t{ a {r)tl a {Q)) and 
their Fourier transforms. We have: 
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F(k,iuj n ) = F\k,iuj n ) = 



Lol + Lot 



and the following notation has been 
(2A fc ) 2 , with A fc = Js 2 /2cosfc 



where v n = 2™L 
introduced: lu 2 = 
and Afc = Jj_/4 — n + Js 2 cos k, recovering the dispersion 
relation predicted by Gopalan, Rice and Sigristc3. As 
shown in Refo, the parameters fj, and s are determined 
by solving the self-consistent saddle point equations. Let 
us now turn to the calculation of the Raman intensity. 

The Raman intensity is proportional to Im^ij(iw„ — > 
lu + iO) where: 

XR^n) = W dre^ (T T (tit a )(r)(tlt p )(0)) (18) 

By using the definition ( ^8[ ) and applying Wick's the- 
orem, the following expression for the Raman suscepti- 
bility is obtained: 

Xfl(iWn) = fl~ l ^2 J ^ [G(k, iv n )G{k, iv n - iu n ) 

+F(k,iv n )F\k,ioj n - iv n )) . (19) 

Performing the usual linear response calculation! 2 !, we 
obtain as a final result: 
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The Raman scattering spectra will show two peaks, one 
at energy u> = 2lu v = 2A S corresponding to the bottom 
of the triplet band, and a second one at lu = 2wo, cor- 
responding to the top of the triplet band. Close to the 
critical frequency lv*, I(lu) ~ (lu — lu*)^ 1 / 2 . This bekavior 
can be easily understood by a density of states argument. 
The resulting spectra is plotted in figure [|. No signal is 
obtained for lu < 2A S in agreement with numericsES. Let 
us note that in recent experiments, a Raman, scattering 
peak at twice the gap is observed in CaV205t] where the 
spin-gap and the exchange constant are estimated to be 
A s ~ 400cm- 1 and Jj_ ~ 640AT. These results are m 
qualitative agreement with our theory. In the case of 
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Sri4 Cu2404ii the situation is more complicated due to 



W. Brenig, Phys. Rev. B 56, 2551 (1997), |cond 



the coexistence in the structure of diipqrizcd spin chains 



having a spin gap A c h a in. = 12meV£2 and-pf spin lad- 
ders having a spin gap Ai a( jdcr. — 32 meVEZI. In Ref. 0, 
a peak was obtained at 570cm -1 ~ 71meV in Raman 
scattering experiments on polycrystalline samples. Ac- 
cordind to our theory, this would lead to a spin gap of 
~ 35meV, in agreement with neutron scattering datas. 
A more recent investigation or Raman scattering on sin- 
gle crystalsu identifies a peak at 498cm -1 as the Raman 
peak associated with the gap. The peak at 569cm -1 is 
identified with a (0, 0) gap. According to the authors of 
Ref. ^, the other peaks are associated with bound states 
or single magnon light scattering. It is known that bound 
states of magnetic excitations can be formed below the 
gap in a spin ladderr 8 tr 9 l. In our treatment, we have 
been neglecting them altogether. They should give rise 
to peaks below_the threshold 2 A, as has been observed 
in experimentsu. This problem is under investigation. 
To summarize, we have considered Raman scattering in 
a spin ladder both in the weak coupling and the strong 
coupling approximation. We have shown that only the 
strong coupling treatment gave rise to peaks in the Ra- 
man intensity. Future directions include the considera- 
tion of the effect of bound states on the Raman spectra. 
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FIG. 1. Raman intensity in arbitrary units for Ji < J at 
T = OK obtained from the Majorana fermion approach. The 
frequency u is measured in units of the gap. 
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FIG. 2. Raman intensity in arbitrary units for J±/J = 2 
at T = OK. The frequency u is measured in units of the gap. 
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